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Abstract 

The polytropic gas model is investigated as an interacting dark energy scenario. 
The cosmological implications of the model including the evolution of EoS parameter 
w\, energy density f^A and deceleration parameter q are investigated. We show that, 
depending on the parameter of model, the interacting polytropic gas can behave as a 
quintessence or phantom dark energy. In this model, the phantom divide is crossed 
from below to up. The evolution of q in the context of polytropic gas dark energy 
model represents the decelerated phase at the early time and accelerated phase later. 
The singularity of this model is also discussed. Eventually, we establish the corre- 
spondence between interacting polytropic gas model with tachyon, K-essence and 
dilaton scalar fields. The potential and the dynamics of these scalar field models are 
reconstructed according to the evolution of interacting polytropic gas. 
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I. INTRODUCTION 



Recent cosmological observations obtained by SNe la WMAP [2|, SDSS [3( and X-ray 
(4] experiments reveal that our universe expands under an accelerated expansion. In the 
framework of standard Freidmann-Robertson- Walker (FRW) cosmology, a missing energy 
component with negative pressure dubbed dark energy (DE) is responsible for this expan- 
sion. The nature of DE is still unknown and scientists believe that the problem of DE is 
a major puzzle of modern cosmology. Up to now, many theoretical models have been in- 
vestigated to interpret the behavior of DE. The time-independent cosmological constant, A, 
with EoS parameter w = — 1 is the earliest and simplest candidate of DE. The cosmological 
constant suffers from two well known difficulties namely "fine-tuning" and "cosmic coinci- 
dence" problems. The alternative candidates for DE problem are the dynamical dark energy 
scenario with time varying EoS parameter ,w. According to some analysis on the SNe la 
observational data, it has been shown that the time-varying DE models give a better fit 
compare with a cosmological constant [5|]. There are two different categories for dynamical 



essence 



[9], tachyon flO|, dilaton H| and so forth, (ii 



DE scenario: (i) The scalar fields including quintessence p|, phantom [7], quintom [8j, K 



a, 



Chaplygin gas models [12j,ll3j, braneworld models [14j |. holographic [15] and agegraphic [1 



The interacting DE models including 



models. The holographic DE model is constructed in the light of holographic principle of 



quantum gravity 



171 | and the agegraphic model is constructed based on the uncertainty re- 



lation of quantum mechanics together with the gravitational effect in general relativity 18 ]. 
The interaction between DE and dark matter is supported by recent observations prepared 
by the Abell Cluster A586 [19J . However the strength of this interaction is not clearly iden- 



tified 



20| | . Also, recent astronomical data supported that our universe is not a perfectly flat 



and has a small positive curvature 21] . 

The polytropic gas model has an important role in stellar astrophysics. It can explain 
the equation of state of de gene rate electrons and degenerate neutrons in white dwarfs and 
neutron stars, respectively [28]. This model can also be useful when the pressure and den- 
sity are adiabatically related to each other in main sequence stars 28]- Here we consider 
the interpretation of dark energy scenario with the EoS parameter of polytropic gas. U. 
Mukhopadhyay and S. Ray used some dynamical A model with polytropic equation of state 
in dark energy scenario 22j. Recently, by using the polytropic gas model, the interaction 
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24J. Karami, et al. obtained the phantom 
Also, karami, et al. reconstructed the 



between DE and dark matter is investigated 
behavior of interacting polytropic gas model 
/(T)-gravity from the polytropic gas DE model 271] . They also studied the correspondence 
between the interacting new agegraphic dark energy model with polytropic gas model in 
non-flat FRW universe and reconstructed the potential and the dynamics for the scalar field 
of the polytropic model to describe the accelerated expansion of the universe {23)]. The 
above statements motivate us to consider more cosmological implications of this model in 
dark energy scenario. One of the interesting features of this model that we discuss is that 
in the polytropic gas dark energy scenario the phantom regime can be achieved even in the 
absence of interaction between dark energy and dark matter. This makes it distinguishable 
from many other DE model whose W\ can not crosses the phantom regime without the 
interaction between DE and dark matter. We consider the interacting polytropic gas as 
a phenomenological DE model. In the phenomenological models of DE the pressure p is 
given as a function of energy density p, i.e., p = —p — f(p) {25]. Considering f(p) = 0, 
the EoS parameter of phenomenological models cross w = —1, i.e., the EoS of cosmological 
constant. Nojiri, et al. investigated four types singularities for some illustrative examples 



of phenomenological models 25(. The polytropic gas model has a type III. singularity in 
which the singularity takes place at a characteristic scale factor a s . 

Here, we obtain the deceleration parameter q to explain the decelerated and accelerated 
expansion phases of the universe dominated by polytropic gas dark energy fluid. The be- 
havior of interacting polytropic gas in the quintessence regime is also calculated. We study 
the correspondence between the tachyon, K-essence and dilaton fields with the interacting 
polytropic gas dark energy and reconstruct the potential and the dynamics of these scalar 
fields according the evolutionary form of interacting polytropic gas model. 

II. POLYTROPIC GAS DE MODEL 

The equation of state (EoS) of polytropic gas is given by 

Pa = Kp\ + \ (1) 



where K and n are the polytropic constant and polytropic index, respectively 28]. 
Assuming a non-flat Friedmann-Robertson- Walker (FRW) universe containing DE and CDM 
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components, the corresponding Friedmann equation is as follows 

where H is the Hubble parameter, M p is the reduced Plank mass and k = 1,0,-1 is a 
curvature parameter corresponding to a closed, flat and open universe, respectively. p m and 
Pa are the energy density of CDM and DE, respectively. Recent observations support a 



closed universe with a tiny positive small curvature Qk =~ 0.02 29]. 
The dimensionless energy densities are defined as 

q _ Prn _ Pm q _ PA _ PA q _ k /q\ 

m ~ p c ~ 3M 2 tf 2 ' p c 3M 2 tf 2 fc ~ a 2 tf 2 W 

Therefore the Friedmann equation (j2J) can be written as 

fi m + A = 1 + ft fc . (4) 

Considering a universe dominated by interacting polytropic gas DE and CDM, the total 
energy density, p = p m + Pa, satisfies a conservation equation 

p + 3H(p + p) =0 (5) 

However, by considering the interaction between DE and dark matter, the energy density of 
DE and dark matter does not conserve separately and in this case the conservation equations 
are given by 

p m + 3Hp m = Q, (6) 
p A + 3H(p A + p A ) = -Q, (7) 

where Q indicates the interaction between DE and CDM. Three forms of Q which have been 
extensively used in the literatures are |3o| 

Q = 3aHp A , 3(3Hp m , 3-fH(p A + p m ), (8) 

where a, ft and 7 are the dimensionless constants. The Hubble parameter H in the Q-terms 
is considered for mathematical simplicity. Indeed, the interaction forms in Eq.Q are given 
by hand, since the Q in Eqs.([6l [7]) should be as a function of H multiplied with energy 
density. Similar to the standard ACDM model, in which the vacuum fluctuations can decay 
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into matter, here the interaction parameter Q indicates the decay rate of the polytropic 
gas into CDM component. Recently, the the interaction between DE and dark matter is 



3l| . For mathematical simplicity, we consider the first form of interaction 



presented in ^ 
parameter Q. 

Using Eq.flTJ, the integration of continuity equation for interacting dark energy component, 
i.e. Eq.((7]), obtains 

where B is the integration constant, K = and a is the scale factor. Note that to have 
a positive energy density for an arbitrary number of n, it is required Ba 3 ^ 1+a ^ n > K. It is 
worthwhile to note that the phantom behavior of interacting polytropic gas has been also 
studied in [24|. In the case of Sffl 3 ' 1 ^/" = K, we have p — > oo and therefore the polytropic 
gas has a finite-time singularity at a c = (K / B) n ^ 3( - 1+a \ This type of singularity, in which at 
a characteristic scale factor a s , the energy density p — > oo and the pressure density \p\ — > oo, 
is indicated by type III singularity [251 ]. 
Substituting Q = 3aHp A in (I7j), we have 

PA + 3#(1 + a + w A )p A = 0, (10) 

Taking the derivative of Eq. (jHJ) with respect to time, one can obtain 

p\ = -WH(l + a)a^pi + " (11) 

Substituting Eq.l flT]) in (flQj) and using Eq.(j9]) , we can obtain the EoS parameter of inter- 
acting polytropic gas as 

3(l + a) 
Ob n 

W A = "1 z^-U (12) 

c — a n 

where c = K/B. By defining the effective EoS parameter as w A = w\ + a = 
— 1 — a ( n ' j(c — a ( n ) ) , we see that the interacting polytropic gas model behaves as a 
phantom model, i.e. w\ < —1, when c > a 3 ( 1+a )/ n . The phantom behavior of polytropic 
gas is similar to generalized chaplygin gas model, where it has been shown that the 
generalized chaplygin gas with negative value of model parameter can behave as a phantom 
dark energy 13( . Note that in the case of phantom polytropic gas, from Eq. ([9]) we see that 
only even numbers of n should be chosen to have a positive energy density. The interesting 
feature of polytropic gas model is that it can obtain the phantom regime even in the 
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absence of interation. For this aim, it is enough to insert a = in Eq. (fl2|) and see that for 
c > a ( n ' the phantom regime, w\ < —1, can be achieved. This makes it distinguishable 
from many other dark energy models whose w\ cannot cross the phantom regime without 
interaction term. The other interesting aspect of the polytropic gas is that the interacting 
polytropic gas dark energy crosses the phantom divide from wa < — 1 to w\ > — 1 (see 
Fig.(l), left panels). This behavior of polytropic gas is similar to interacting agegraphic 
dark energy model in which the phantom divide is crossed from below to up (see Figs. (2,3) 
of [261 )• The similarity of the interacting agegraphic dark energy and polytropic gas is that 
both models cross the phantom divide from below to up. 

The interacting polytropic gas behaves as a quintessence model, i.e. —1 < w e / < —1/3, 
when -oo < c < -a 3(1+a)/n /2. The condition -a 3(1+a)/n /2 < c < a 3 ( 1 + Q )/™ leads to 
ufjl* > — 1/3 and consequently the accelerated expansion, in this case, can not be achieved. 
At c = a , 3 ( 1 + a )/™ ; the interacting polytropic gas has a singularity. Hence, depending on 
the parameter c, the polytropic gas can behaves as a phantom or quintessence models of 
DE. Also it is worth to mention that the polytropic gas model behaves as a cosmological 
constant, i.e.,u^ — > —1, at the early time (i.e. a — > 0) whereas the universe is dominated 
by pressureless dark matter. 

In Fig. (1), the evolution of w\ as a function of scale factor is plotted for different values 
of the parameters c and n.The other interesting aspect of the polytropic gas is that the 
interacting polytropic gas dark energy crosses the phantom divide from w\ < —1 to 
wa > —1 (see Fig.(l), left panels). This behavior of polytropic gas is similar to interacting 
agegraphic dark energy model in which the phantom divide is crossed from below to up (see 
Figs. (2,3) of 26) ). The similarity of the interacting agegraphic dark energy and polytropic 
gas is that both models cross the phantom divide from below to up. In upper panels we 
fix the polytropic index as n = 2 and in lower panels the parameter c is fixed. In upper 
left panel the negative values of c are selected to obtain the transition from phantom to 
quintessence regime. In upper right panel, the positive values of c are selected. In this case 
the interacting polytropic gas behaves as a phantom like field. Same as left panel, we fix 
the polytropic index n — 2. Here, one can easily find the phantom behavior of polytropic 
gas model. It is worth noting that the phantom regime of polytropic gas model is restricted 
with a characteristic scale factor a s = c n / 3 ( 1 + a ) ; where we encounter with a singularity 
at this epoch. In lower panels of Fig.(l), the dependency of the evolution of wa on the 
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polytropic index parameter n is studied. In lower left panel, by fixing c = —1, we studied 
this dependency for polytropic gas model. It is easy to see that the larger value of n gets 
the larger w\ at a < 1 and smaller w\ at a > 1. In lower right panel, by fixing c = 2, the 
dependency of wa on the parameter n is investigated for phantom polytropic gas model. 
Unlike to lower left panel, the larger value of n gets the smaller w\ at a < 1 and larger wa 
at a > 1. 



In order to obtain the evolution of dimensionless energy density, ^a, let us start with 
Eqs.Q and ([3]) and obtain the density parameter of interacting polytropic gas as 

3(l + a) ~ 

Qa ~ SMjlP (13) 

Taking the derivative of Eq.f fTBl with respect to time and using Q' — Q/H, we can obtain 

K = M^i^~^ +2 w) < 14 > 

v a « — c n 7 

where prime denotes the derivative with respect to x — In a. Taking the derivative of 
Friedmann equation ([2]) with respect to time and using Eqs.([9]), (B|, (J6]), ( |T3l) and Q = 
3aHpk, one can find that 



H_ 3 

IP ~ ~2 



c(l + a) Q k 

"A ~3(i+^ + 1 + "5- 

a n — c 



(15) 



Substituting this relation into Eq. (TH]) . we obtain the evolutionary equation for energy den- 
sity parameter of interacting polytropic gas as: 



a 3(1 + a) _ C Q Q 

tt' A = -3tt A - ^^y (1 - Ov) + a — 3^ - -y , (16) 



a n — c a n — c 

where flk is given by 

tt k = a>y- — - (17) 

1 — 07 

and 7 = Q k0 /Q m0 . 

In Fig. (2), by solving the differential equation ( fl6l) . we show the evolution of fl\ for different 
model parameters c and n as well as different interaction parameter a. Here we assume only 
the positive values of c, i.e., the phantom polytropic gas model. The numerical values of 
density parameters at the present time are taken as: ^ao — 0.7, Q m Q = 0.3 and Q k o — 0.02. 
In upper panels, we consider the non-interacting polytropic gas and in lower panel the 
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interaction term is included. Here, we see that £l\ — > at the early time and tends to 1 at 
the late time. Hence the polytropic gas model can describe the matter-dominated universe 
in the far past. Also, at the late time, we encounter with dark energy dominated universe 
(Qa — > 1). In upper left panel, by fixing the parameter n, the polytropic gas starts to be 
effective earlier and tends to a lower value at the late time when c is larger. On the other 
hand, in upper right panel, we see that for fixed parameter c, the polytropic gas starts to 
be effective earlier and also tends to a higher value at the late time when n is smaller. 
In lower panel, the effect of interaction parameter a on the evolution of ^a is studied. Here 
one can see that the polytropic gas starts to be effective earlier, by increasing the interaction 
parameter a. Also, at a > 1, the parameter Qa is smaller for larger values of a. 
For completeness, we derive the deceleration parameter 



for polytropic gas model. Substituting Eq.( fl5|) in ( fl8l) we get 



It is worth noting that in the limiting case of matter-dominated phase and considering flat 
universe in the absence of interaction term, Eq. flT9|) is reduced to q — 1/2 which represents 
the decelerated expansion (q > 0) of the universe. 

In Fig. (3), we show the evolution of q as a function of a for different model parameters 
c and n as well as different interaction parameter a. Here we discuss the evolution of q 
for phantom polytropic gas model, by assuming positive c. Upper panels is plotted in the 
absence of interaction between dark energy and dark matter and the lower panel is plotted 
in the presence of interaction term. The parameter q converges to 1/2 at the early time, 
whereas the universe is dominated by pressureless dark matter. In upper left panel, by 
fixing n, the accelerated expansion is achieved earlier by increasing c. Also, in the upper 
right panel, we see that by increasing n, q becomes larger at the deceleration phase and 
gets smaller at the acceleration phase. It is worth noting that, although, both the model 
parameters n and c impact the evolution of deceleration parameter q, but the change of sign 
from q > to q < depends on the parameter c of the model. 




a 



— c 



(19) 
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III. CORRESPONDENCE BETWEEN POLYTROPIC GAS DE MODEL AND 

SCALAR FIELDS 

In the present section we establish a correspondence between the interacting polytropic 
gas model with the tachyon, K-essence and dilaton scalar field models. The importance 
of this correspondence is that the scalar field models are an effective description of an 
underlying theory of dark energy and therefore it is worthwhile to reconstruct the potential 
and the dynamics of scalar fields according the evolutionary form of polytropic gas model. 
For this aim, first we compare the energy density of polytropic gas model (i.e. EqJH]) with 
the energy density of corresponding scalar field model. Then, we equate the equations of 
state of scalar field models with the EoS parameter of polytropic gas (i.e. EqJT2l). 



A. Polytropic gas tachyon model 

It is believed that the tachyon can be assumed as a source of DE [33]. The tachyon is an 
unstable field which can be used in string theory through its role in the Dirac-Born-Infeld 



(DBI) action to describe the D-bran action 34|. The effective Lagrangian for the tachyon 
field is given by 

where V(<f>) is the potential of tachyon. The energy density and pressure of tachyon field 



arc 



34] 



F =, ,20) 



1 -(f) 2 

p*=-v(0)\A-^ 2 - ( 21 ) 

The EoS parameter of tachyon can be obtained as 

«V = ^ = 2 -1. (22) 
P4> 

In order to have a real energy density for tachyon field, it is required that — 1 < <p < 1. 
Consequently, from Eq. (122ft . the EoS parameter of tachyon is constrained to — 1 < < 0. 
Hence, the tachyon field can interpret the accelerates expansion of universe, but it can not 
enter the phantom regime, i.e. w\ < — 1. In order to reconstruct the potential and the 
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dynamics of tachyon according to evolution of interacting polytropic gas model, we should 
equate Eqs.flHD and © and also Eq.® with Eq.fl20D as follows 

3(l+c) 

(X n • n 

w A = -1 3(T+^r - « = - 1- (23) 

c — a « 

*-(j^-t)"-T^ (-) 



Hence we get the following expressions for dynamics and potential of tachyon field 

3(l+a) 

a (25) 



3(l+a) 

c — a n 



( 1 V 

= \l 1 + s^M + a - 3 (1+a) ~ (26) 

c — a n \Ba « — K / 

For c > a 3( - 1+a - )//n , from Eq. (l25|) . we obtain 2 < which represents the phantom behavior 

of tachyon field. It is worth noting that the reconstructed tachyon field according to the 

interacting polytropic gas can cross the phantom divide. By definition <p = iip and changing 

the time derivative to the derivative with respect to logarithmic scale factor, i.e. d/dt = 

Hd/dx, the scalar field ip can be integrated from Eq.( |25l) as follows 



3(l+a) 

■ CL 71 

V>(a:) - <0(O) = / 77 V MT^)~ adx (27) 

o n V c — a n 

B. Polytropic gas K-essence model 

The idea of the K-essence scalar field was motivated from the Born-Infeld action of string 



theory and can explain the late time acceleration of the universe 35| . The general scalar 
field action for K-essence model as a function of <p and \ = 2 /2 is given by [36( 

S = J d'x^-g p(<j>,x), (28) 

where the Lagrangian density p(<p, x) relates to a pressure density and energy density through 
the following equations: 

p(<P,x) = fm-X + X 2 ), (29) 
p(0,x) = /(0)(-X + 3x 2 )- (30) 
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Hence, the EoS parameter of K-essence scalar field is obtained as 



By comparing Eqs.(IT2l and ( l3il . we have 



1 



3(l+a) 

a n x 
— o = 

3(l + a) 

c — a n 



3X 



(31) 



(32) 



Hence the parameter x is obtained as 



X 



3(l+a) 

2 + 3("l+c) + OL 

c—a n 
3(l+a) 

4 + 3 g 3( " 1+a) +3a 

c—a n 



(33) 



From Eq. (13 II) . one can see the phantom behavior of K-essence scalar field (wk < — 1) when 
the parameter x hes in the interval 1/3<x<1/2. 

Using <fi 2 = 2x and changing the time derivative to the derivative with respect to x — In a, 



we obtain 



1 

H 



\ 



4 + 2- 



3(l+a) 
I n 
3(l+c) 



+ 2a 



c—a n 



3(l + a) 

4 + 3 a 3( " 1+a) +3a 



(34) 



c—a n 



The integration of Eq. (j3^|) yields 



0(x) - 0(0) 



i 



3(l + a) 

4 + 2 g 3( " 1+a) +2a 

c—a « 



4 + 3^ 



3(l + a) 



(35) 



3(l + a) 



+ 3a 



Here, we reconstructed the potential and the dynamics of K-essence scalar field according to 
the evolutionary form of the interacting polytropic gas model. The K-essence polytropic gas 
model can explain the accelerating universe and also behaves as a phantom model provided 
1/3 < X < 1/2. 



C. Polytropic gas dilaton model 

A dilaton scalar field can also be assumed as a source of DE. This scalar field is originated 
from the lower-energy limit of string theory (3^ . The dilaton filed is described by the effective 
Lagrangian density as 

PD = -X + ce x<t> x 2 , (36) 
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where c and A are positive constant. Considering the dilaton field as a source of the energy- 
momentum tensor in Einstein equations, one can find that the Lagrangian density corre- 
sponds to the pressure of the scalar field and the energy density of dilaton field is also 
obtained as 

PD = -x + 3ce A V, (37) 

Here 2\ = <fi 2 - The negative coefficient of the kinematic term of the dilaton field in Einstein 
frame makes a phantom like behavior for dilaton field. The EoS parameter of dilaton is 
given by 

Pd -l + ce x ^x 
UD -p^- -l + 3ce»x { ) 

In order to consider the dilaton field as a description of polytropic gas, we establish the corre- 
spondence between the dilaton EoS parameter, wd, and the EoS parameter wa of polytropic 
gas model. By equating Eq. fl38l) with Eq. ffT2"]) . we find 

3(l + e») 

2 + a ^ + a 



w A -l " 1 £ 

4 + 3^^— + :{n 



ce ^ x = WA = £z^L (39) 

A 3w A - 1 



3(l+c) 
c—a ri 



By using \ = <f) 2 /2 and <fi — 4>'H, the scalar field (f) can be obtained as 



2 ln (>(0)/2 + A 



A V V2c 7n H 



\ 



3(l+ct) 

2 + 3 n + „i + a 



c—a n 



3(l+a) 

dx) (40) 



3(l + a) 

4 + 3 a 3(" +a) +3a 

c—a n 



Here we presented the reconstructed potential and dynamics of dilaton scalar field according 
to the evolution of interacting polytropic gas model. 



IV. CONCLUSION 



In this work we presented the interacting polytropic gas model of dark energy to 
interpret the accelerated expansion of the universe. Assuming a non-flat FRW universe 
dominated by interacting polytropic gas DE and CDM, we studied the cosmic behavior of 
polytropic gas model. For this aim, we calculated the evolution of effective EoS parameter 
and showed that for positive values c > a 3 ( 1+Q )/ ra with even numbers of n this model behaves 
as a phantom DE model and in the case of — oo < c < — a 3 / n /2 it treats as a quintessence 
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model. Similar to interacting agegraphic dark energy model, the interacting polytropic gas 
model crosses the phantom divide from below (w\ < —1) to up (w\ > — 1). The transition 
from phantom to quintessence depends on the parameter c of the model. For larger value 
of c, the transition take place sooner. In the case of phantom polytropic gas model, w\ is 
larger for larger value of c. In the scenario of polytropic gas model, the phantom divide can 
be crossed even in the absence of interaction. We also calculated the evolution of energy 
density f2 A . The matter dominated phase at the early time and DE-dominated universe at 
the late time can be described in the context of polytropic gas model. The polytropic gas 
starts to be effective earlier for larger value of interaction parameter as well as for larger 
value of the parameter c or smaller value of n. We calculated the deceleration parameter 
q and obtained the decelerated and accelerated expansion phases of the universe in the 
context of polytropic gas model. The transition from decelerated expansion (q > 0) to 
accelerated expansion (q < 0) takes place sooner for larger value of c and also by increasing 
the interaction parameter a. Since the scalar fields models are the underlying theory of 
dark energy, we proposed a correspondence between interacting polytropic gas model with 
the tachyon, K-essence and dilaton scalar fields models. We reconstructed the potential and 
the dynamics of these scalar fields according to the evolution of the interacting polytropic 
gas model. 
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FIG. 1: The EoS parameter w/\_ of interacting polytropic gas model as a function of cosmic scale 
factor a for different model parameters c and n. The interaction parameter is chosen as a = 0.1. 
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FIG. 2: The evolution of energy density ,0a, in terms of cosmic scale factor a for interacting 
polytropic gas model. In upper left panel, the parameter n is fixed and the parameter c is varied. 
In upper right panel, we fix c and vary n. In lower panel, by fixing the parameters c and n, we 
vary the interaction parameter a. 
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FIG. 3: The evolution of deceleration parameter q versus of scale factor a for interacting polytropic 
gas model. In upper left panel, the parameter n is fixed and the parameter c is varied. In upper 
right panel, we fix c and vary n. In lower panel, by fixing the parameters c and n, we vary the 
interaction parameter a. 
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